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Abstract 



We present a computation of the 0(a s ) QCD corrections to W Zj production 
at the Large Hadron Collider. The photon is considered as real, and we include full 
leptonic decays for the W and Z bosons. Based on the structure of the VBFNLO 
program package, we obtain numerical results through a fully flexible Monte Carlo 
program, which allows to implement general cuts and distributions of the final-state 
particles. The NLO QCD corrections are sizable and strongly exceed the theory error 
obtained by a scale variation of the leading-order result. Also, the shapes of relevant 
observables are significantly altered. 
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Figure 1: Examples of the three topologies of Feynman diagrams contributing to the process 
pp -^■W ± Z r ) J r X at tree-level. In the right-hand diagram, the quartic coupling is marked 
with a dot. 

1 Introduction 

At the LHC a new energy frontier is reached, which allows for new searches of unknown 
particles and further tests of the Standard Model. From the theory side this also requires 
very precise predictions of cross sections and distributions, exceeding the precision of a 
leading-order approximation. We present a calculation of W ± Zj production including QCD 
corrections and the leptonic decays of the W- and the Z-boson with off-shell effects, 



This process with leptons, photons and missing energy in the final state provides a back- 
ground to new physics searches (see, for example, Ref. pQ). Also, this process offers the 
possibility to study the quartic vector-boson couplings WW Z^ and WWjj (see right dia- 
gram in Fig. [1]) [2] and test the Standard Model. It is one of the missing pieces for a full 
knowledge of triple vector boson production at next-to leading order (NLO) QCD. There 
has been a strong effort for the calculation of these processes. The processes with only mas- 
sive vector bosons in the final state have been completely calculated [3HB] . Rather recently, 
also the NLO QCD calculation of the processes WWj and ZZj have been completed [7] 
and first results of the NLO QCD computation of the process W'yy have been presented [8]. 
Vector-boson pair production accompanied with one jet has also been studied including QCD 
corrections for WWj, W73, WZj and ZZj [9HI3]. 

2 Calculational Details 

We calculate all contributions to the processes (11.11) up to order a s a 5 in the limit where 
all fermions are massless. At leading order, 71 distinct diagrams appear, which we group as 
three different topologies, according to the number of vector bosons attached to the quark 
line. An example of each class is depicted in FigJT] For the last two topologies we also 
include the cases where vector bosons are radiated off the lepton lines. 

To speed up the calculation, invariant subparts, which appear multiple times in different 
Feynman graphs, are computed only once per phase-space point and independently of the 



pp,pp -»■ W ± Z 1 + X -»■ if W i+ ^2 7 + X . 



(1.1) 
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rest of the cross-section. Hereby, we use the procedure of leptonic tensors, as first described 
in Ref. [14] . This greatly reduces the computation time needed. For the computation of 
the matrix elements, we use the helicity method introduced in Ref. fTS]. Furthermore, by 
charge conservation, a W boson must always couple to the quark line. Hence, we only need 
to compute the left-handed chirality part. 

At NLO QCD, virtual and real emission diagrams contribute to the cross section. Both 
contain infrared divergences, which must cancel in the sum according to the Kinoshita-Lee- 
Nauenberg (KLN) theorem [16] . To handle this cancellation in a numerically stable way, we 
use the Catani-Seymour dipole subtraction method [T7]. Initial-state collinear singularities 
are partly factorized into the parton-distribution functions. This leads to additional so-called 
"finite collinear terms". 

The NLO real corrections are given by diagrams where an additional gluon is attached 
to the quark line which is possible in two different ways. Either this gluon is a final-state 
particle and considered as radiated off the quark line, or an initial-state gluon which splits 
into a qq pair and we have an emission of a quark. With 194 different Feynman diagrams, 
the use of leptonic tensors proves to be an advantage in this case. 

The presence of isolated on-shell photons requires extra care in the case of real emissions. 
An additional singularity arises from photon emission collinear to a massless quark. Requir- 
ing a simple separation cut between photon and jet is not allowed, since the cancellation 
of the gluonic infrared divergences between virtual and real emission processes would be 
spoiled. The phase space for soft-gluon emission would be reduced while leaving untouched 
the virtual counterpart, leading to the non-cancellation. In principle, this problem can be 
solved by adding processes with quark fragmentation to photons. Here, we use a simpler 
approach by requiring a specially crafted cut first suggested in Ref. [18]. We will discuss the 
details of this cut in section |3j 

Virtual NLO QCD corrections arise from the insertion of gluon lines into the topologies 
of Fig. [T] in every possible way. Therefore, the left topology gives rise to loops with up to five 
external particles, i.e. pentagon diagrams. For the middle and right ones at most box and 
triangle diagrams, respectively, appear. Up to the box level, we compute the loop integrals 
using Passarino-Veltman reduction [19] ; additionally we avoid the explicit calculation of the 
inverse of the Gram matrix. Instead, we solve a system of linear equations, which is numer- 
ically more stable close to the singular points. For the pentagon contributions, we apply the 
method of Ref. [20]. This circumvents the appearance of small Gram determinants in planar 
configurations of the external momenta altogether. The complete virtual corrections 

M B , (2.1) 

can be separated into a part which is proportional to the Born matrix element, Mb, and 
a remainder, My. Here, Q denotes the partonic center-of-mass energy, which corresponds 
to the invariant mass of the WZ'y final state. For diagrams with only a single W attached 
to the quark line, My vanishes, i.e. the virtual corrections completely factorize to the Born 
amplitude. 
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For this factorization formula to hold, the transversality property of the photon e(p 7 )-p 7 = 
must be used j7]. Then we can reuse the general results for the finite remainder already 
obtained in Refs. [SIE], adapting the attached vector bosons to our case. Another opti- 
mization can be performed by shifting parts from the pentagon diagrams into box contribu- 
tions I1JEJC0EE1] • To do this, we split the polarization vector of the vector bosons into a part 
proportional to the four-momentum of the boson qy and a remainder e v , which is chosen 
such that 

iv ■ (qw + qz) = 0. (2.2) 

This part leads to a reduction in the size of the pentagon contributions, which we can 
therefore compute with lower statistics. Contracting the pentagons with q v , they can be 
expressed as a difference of two box diagrams which are numerically faster to compute. 
Hence, we can gain speed while keeping the total error the same. This shift also serves as 
a consistency check between the box and pentagon routines, as the total result must stay 
unchanged. 

To verify the correctness of our calculation, we have performed several checks. First, we 
have compared all tree-level amplitudes against matrix elements generated by MadGraph [2T] 
and find an agreement of at least 14 digits, which is at the level of the machine precision. 
Additionally, we have also compared the integrated cross section for both pp — >■ W ± Z^ and 
pp — > W ± Z'-fj against MadEvent and Sherpa [22]. We find an agreement at the per mill 
level, which is compatible with the integration error. Furthermore, we have checked the 
implementation of the Catani-Seymour subtraction scheme. We have verified that for the 
real emission part the ratio between the differential real-emission cross section and the dipole 
approaches —1 once we go to the soft or collinear limit. Also, we have checked that finite 
contributions left after the cancellation of the infrared divergence can be shifted between the 
virtual and real parts without affecting the total result. 

For the phase space implementation, the 2—7-5 phase space was split into three different 
parts. They correspond to three topologies which we have identified to give the largest 
contributions, which are of comparable size each. The first one corresponds to a 2 — > 2 s- 
channel structure, where possibly a massless particle (the final-state jet for the real-emission 
contribution) is radiated off first. The two outgoing particles are the massless photon and a 
pseudo-particle with a mass of 450 GeV and a width of 100 GeV. The latter one then decays 
into the two massive vector bosons, which in turn decay into their respective fermion-anti- 
fermion pair, all via 1 — > 2 topologies. For all invariant masses of intermediate particles 
we apply a Breit-Wigner mapping to the corresponding random number, so that the peak 
is flattened out, integrating over the whole available energy range. The second and third 
structure correspond to 2 — > 2 s-channel production of the two massive vector bosons, where 
one undergoes a two-body decay into its fermions and the other one, via a 1 — > 3 process, into 
the fermions and the photon. To avoid double-counting we compute for both fermion-anti- 
fermion pairs the difference of the invariant mass of the photon and the fermion-anti-fermion 
pair and the mass of the corresponding vector boson. If both differences are larger than 30 
GeV, the phase-space point is assigned to the first topology, otherwise to that one where the 
difference is smaller. 
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3 Results 



We perform the numerical evaluation of our calculation with an NLO Monte Carlo pro- 
gram based on the structure of the VBFNLO program package [23]. As input parameters in 
the electroweak sector we take the W and Z boson masses and the Fermi constant. The 
weak mixing angle and the electromagnetic coupling constant are computed from these using 
tree-level relations: 

m w = 80.398 GeV sin 2 (6 W ) = 0.22264 

m z = 91.1876 GeV a" 1 = 132.3407 

G F = 1.16637 • 10~ 5 GeV" 2 . (3.1) 



Top-quark effects are not considered and all other quarks are taken massless. Effects 
from generation mixing are neglected, as we set the CKM matrix to the identity matrix. As 
the central value for factorization and renormalization scales we choose the invariant mass 
of the leptons and the photon 



fJ-F = = A*o = rn WZj = yj (p<! + Vvx + Ph + Ph + Pi) 2 - ( 3 - 2 ) 

For the parton distribution functions, we choose CTEQ6L1 at LO and the CTEQ6M set 
with a s {m z ) = 0.1176 at NLO [21]. 

We impose the following set of minimal cuts on the rapidity y and the transverse mo- 
mentum, p T , of the final-state photon and charged leptons 

p Ti > 10 GeV p Te > 20 GeV |y 7 | < 2.5 \y e \ < 2.5 . (3.3) 

These take into account typical requirements of the experimental detectors. Furthermore, 
leptons, photon and jet need to be well separated in order to avoid divergences from collinear 
photons and to be able to identify them as separate objects in the detectors. Therefore, we 
impose 

R ee > 0.3 R h > 0.4 R je > 0.4 R h > 0.4 m u > 15 GeV . (3.4) 

Here, a jet refers to a final-state quark or gluon in the NLO real emission contribution with 
Pt > 20 GeV and \yj\ < 4.5. The last cut in Eq. ( 13 .4)) eliminates the singularity from a 
virtual photon 7* — y by requiring that the invariant mass of each pair of oppositely 
charged leptons is larger than 15 GeV. For treating the collinear singularity between the 
photon and a parton i, we use the procedure of Ref. [18J. The event is accepted only if 

1 — cosi? 7i 

Pt, <Pt~,- r~ 01 R ii > °o ' ( 3 - 5 ) 

1 — cos do 

where 5o is a fixed separation parameter which we set to 0.7. Eq. (I3.5P allows final-state 
partons arbitrarily close to the photon axis as long as they are soft enough. Thereby, it 
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LHC 


y/s =7 TeV 
LO[fb] NLO [fb] K-factor 


v/i =14 TeV 
LO[fb] NLO [fb] K-factor 


Aa^W+Zi 1 -> e+z/ e /i+/i-7) 
p Tj(e) > 10(20) GeV 
pr, K) > 20(20) GeV 


0.650 1.075 1.65 
0.278 0.482 1.74 


1.324 2.441 1.84 
0.587 1.169 1.99 


4(7("W-Z 7 " ->• e-z/ e /i + /i- 7 ) 
p T , m > 10(20) GeV 
Pr,^ > 20(20) GeV 


0.352 0.621 1.76 
0.146 0.275 1.88 


0.886 1.717 1.94 
0.381 0.813 2.13 



Table 1: Total cross sections at the LHC with center-of-mass energies of 7 and 14 TeV for 
pp — > W ± Z^f+X with leptonic decays, at LO and NLO, and for two sets of cuts. Relative 
statistical errors of the Monte Carlo are at the per mill level. The factor 4 accounts for all 
combinations of final-state first- and second-generation leptons. 



Tevatron (^7s = 1.96 TeV) 


LO[ab] NLO [ab] K-factor 


Aa( ,, W ± Z 1 n ->■ e ± u\ /i+/i~7) 
p Tj(e) > 10(10) GeV 
p Ty(i) > 20(10) GeV 


250.0 370.8 1.48 
107.4 160.1 1.49 



Table 2: Total cross sections at the Tevatron for pp — > W + Z-y+X or pp — > W'Z'j+X 
including leptonic decays, at LO and NLO, and for two sets of cuts. Relative statistical errors 
of the Monte Carlo are below the per mill level. The factor 4 accounts for all combinations 
of final-state first- and second-generation leptons. 
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Figure 2: Left: Scale dependence of the total LHC cross section at y/s = 14 TeV for 
pp ->■ W+Z-f + X -> ^4^2 7 + ?4 T + X at LO and iVLO witiin tie cuts of Eqs. J23 ESI 
13.5)) . Tiie factorization and renormalization scales are together or independently varied in 
the range from 0.1 • u. Q to 10 • /Iq. Right: Same as in the left panel but for the different NLO 
contributions at u, F = = £/i with /i = m WZl . 



retains the full QCD pole, which cancels against the virtual part, but it does not introduce 
divergences from the electroweak sector. 

In Tables [Hand El we present results for the integrated cross section of W ± Z r ) production 
for the LHC with a center-of mass energy of both 7 and 14 TeV and for the Tevatron with 
its center-of-mass energy of 1.96 TeV. Note, for the Tevatron, the cross section for W + Z^j 
and W~Z / ~f production is the same; the given number is the individual result of one of them. 
Besides the standard cut on p^ of 10 GeV we also show results for 20 GeV. For the Tevatron, 
we have reduced the cut on p Tl to 10 GeV throughout. The cross sections shown correspond 
to the production of both electrons and muons for all leptons. Interference effects from 
identical leptons in the final state are neglected, since their contribution is small. 

From hereon, we will focus on the LHC with a center-of-mass energy of 14 TeV. In 
Figs. [2] and |3] we show the dependence of the cross section for W + Z / y and W~ Zj production, 
respectively, when varying the renormalization and factorization scale in the interval 

Hf,Hr = £-IM> (0.1<£<10) (3.6) 

around the central scale /iq given in Eq. (13. 2\i . We see that the variation of the LO cross 
section with the scale strongly underestimates the size of the NLO contributions. At the 
central scale, we obtain a K-factor of 1.84 for W + Z r j and 1.94 for W'Zj. The dependence 
on the factorization scale slightly reduces when we move from a LO calculation to NLO as 
expected. On the other hand, the dependence on the renormalization scale shows a large 
variation. This is due to the fact that a s enters the cross section only at NLO, where we 
observe the typical leading renormalization scale dependence. When varying the factorization 
and the renormalization scale jointly by a factor 2 around the central scale /xq, we see a change 
of 7.5% at LO and 6.7% at NLO for W + Z7. For W'Z'j, the numbers change only slightly 
to 7.7% and 7.3% for LO and NLO, respectively. 
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Figure 3: Left: Scale dependence of the total LHC cross section at y/s = 14 TeV for 
pp -> W-Z-f + X -> ^^2"7 + + x + x at L0 and NL0 within the cuts of Eqs. (E3 
\3.4\ 13.5)) . Tie factorization and renormalization scales are together or independently varied 
in the range from 0.1 • /in. to 10 • /in. Right: Same as in the left panel but for the different 
NLO contributions at fip = /iR = £/i with /to = m^ r 



On the right-hand side of Figs. [2] and [31 we show the combined factorization and renor- 
malization scale dependence of the NLO cross section split into the individual contributions. 
Almost the entire scale dependence is given by the real-emission part, which contains the true 
real-emission cross section, the dipole terms from the Catani-Seymour subtraction scheme 
and the finite collinear terms. We obtain the bulk of the NLO contribution from the Born 
matrix element and the virtual corrections proportional to it. This includes the terms from 
boxes and pentagons factored out in Eq. (12. ip . At the central scale, it is more than twice as 
large as the real part. The finite virtual remainders due to box and pentagon corrections, 
which are shifted using Eq. (12. 2p . only yield a small contribution. The shape of their scale 
dependence is similar to the one of the total cross section. 

In Figs. HI |5] and |6j we show the distribution of the transverse momentum of the photon 
and the hardest lepton as well as of the missing transverse momentum originating from the 
neutrino, respectively. On the left-hand side of each figure, we depict the differential cross 
section at LO and NLO both for W + Z , y and W'Z'y, and on the right-hand side we plot the 
differential K-factor, defined in the following way: 

We present results which do not include any cut on the additional jet, as well as including 
a veto on jets with px > 50 GeV. As previously, a jet is defined as a final-state parton with 
|%-|<4.5. 

We see that for each of the distributions the K-factor for the differential cross sections 
without jet veto is not constant, but shows a strong dependence on the momentum scale. In 
all cases it is close to the integrated one for small values of the transverse momenta coinciding 
with the bulk of the cross section. For large transverse momenta we observe much bigger 
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Figure 4: Left: Transverse-momentum distribution of the photon in W + Z^ and W~Z~f 
production with leptonic decays of the W- and the Z-boson at the LHC with y/s = 14 TeV. 
LO (dashed blue line) and NLO (solid red) results are shown for fi F = u, R = /i = m wz -y and 
the cuts of Eqs. / 13.31 \3.4\ \3.5\) . Right: K- factor for the transverse-momentum distribution 
of the photon as defined in Eq. ( 13. 7\) without (solid red) and including a jet veto of 50 GeV 
(dashed blue). 




Figure 5: Left: Transverse-momentum distribution of the lepton with largest transverse 
momentum in W + Z / y and W~Z / y production with leptonic decays of the W- and the Z- 
boson at the LHC with a/s = 14 TeV. LO (dashed blue line) and NLO (solid red) results are 
shown for /ip = fi R = fi = m W z-y and the cuts of Eqs. / I3.3[ 13.41 13.5)) . Right: K-factor for 
the maximal-transverse- momentum distribution of the lepton as defined in Eq. ( 13.7]) without 
(solid red) and including a jet veto of 50 GeV (dashed blue). 
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Figure 6: Left: Distribution of missing transverse momentum in W + Z r y and W~Z^ pro- 
duction with leptonic decays of the W- and the Z-boson at the LHC with yfs = 14 TeV. 
LO (dashed blue line) and NLO (solid red) results are shown for n F — //r = /i — m wz-y 
and the cuts of Eqs. I\3.3\ \3.4[ \3.5\) . Right: K- factor for the missing transverse-momentum 
distribution as defined in Eq. A3.7\) without (solid red) and including a jet veto of 50 GeV 
(dashed blue). 
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Figure 7: Left: Separation of the photon and the hardest lepton in W + Z / y and W~Z~f 
production with leptonic decays of the W- and the Z-boson at the LHC with yfs = 14 TeV. 
LO (dashed blue line) and NLO (solid red) results are shown for fip = fiR = /i = mwz-y 
and the cuts of Eqs. h3.3[ \3.4[ \3.5\) . Right: K- factor for the distribution of the photon-lepton 
separation as defined in Eq. \3.7\) without (solid red) and including a jet veto of 50 GeV 
(dashed blue). 
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K-factors, typically extending up to a value of three. Once we include the additional jet 
veto, this strong dependence is largely removed. For the transverse-momentum distribution 
of the hardest lepton in W~Z r y for example (Fig. [5]), we obtain K-factors between 1.10 
and 1.33 over the plotted momentum range. The large differential K-factors are therefore 
caused by the emission of the additional jet, where the leptonic system recoils against the jet. 
The integrated K-factors are also reduced, namely to 1.23 for W + Z'-f and 1.29 for W~ Z^j 
production. 

In Fig. [7J we show the separation in the rapidity-azimuthal-angle plane (R separation) 
between the photon and the lepton with the largest transverse momentum. Again, we observe 
a significant dependence of the K-factor on the value of the R separation, varying between 
1.4 and 2.55. Once we include the jet veto, this dependence is again much smaller. 

Therefore, a simple approximation of rescaling the leading-order cross section with the 
integrated K-factor does not hold. It is necessary to perform a full NLO calculation to 
reproduce the correct shape of the distributions. 

4 Conclusions 

We have calculated the NLO QCD corrections to the processes pp,pp — > W ± Z r y + X 
including full leptonic decays of the W and Z boson. With three leptons, a photon and 
missing transverse energy as signature, it is an important background for searches for new 
physics, in particular supersymmetry. Additionally, it can serve as a signal process for 
measuring the quartic gauge couplings WW Z^ and WW'yy at the LHC. 

We find that the corrections yield a sizable increase of the cross section with respect to the 
leading-order result, with integrated K-factors typically around 1.9. The LO scale variation 
strongly underestimates these contributions with a variation below the 10% level. Varying 
factorization and renormalization scale fip — = fi by a factor 2 around the central value 
fio = rniyz-y, we obtain a remaining scale dependence at NLO of about 7%. 

The NLO corrections also show a significant dependence on the observable and on dif- 
ferent phase-space regions. Therefore, it is important to have a dedicated fully-exclusive 
NLO parton-level Monte Carlo code available for W ± Z'y production. This process will be 
included into a future version of the VBFNLO program package. 
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The Feynman diagrams in this paper were drawn using Jaxodraw [251 
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